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ABSTRACT 

The long-term (over more than one decade) X-ray emission from two massive stel¬ 
lar systems known to be particle accelerators is investigated using XMM-Newton. 
Their X-ray properties are interpreted taking into account recent information about 
their multiplicity and orbital parameters. The two targets, HD 168112 and HD 167971 
appear to be overluminous in X-rays, lending additional support to the idea that a 
significant contribution of the X-ray emission comes from colliding-wind regions. The 
variability of the X-ray flux from HD 168112 is interpreted in terms of varying sepa¬ 
ration expected to follow the 1/D rule for adiabatic shocked winds. For HD 167971, 
marginal decrease of the X-ray flux in September 2002 could tentatively be explained 
by a partial wind eclipse in the close pair. No long-term variability could be demon¬ 
strated despite the significant difference of separation between 2002 and 2014. This 
suggests the colliding-wind region in the wide orbit does not contribute a lot to the to¬ 
tal X-ray emission, with a main contribution coming from the radiative shocked winds 
in the eclipsing pair. The later result provides evidence that shocks in a colliding- 
wind region may be efficient particle accelerators even in the absence of bright X-ray 
emission, suggesting particle acceleration may operate in a wide range of conditions. 
Finally, in hierarchical triple 0-type systems, thermal X-rays do not necessarily con¬ 
stitute an efficient tracer to detect the wind-wind interaction in the long period orbit. 

Key words: Stars: early-type - Stars: binaries - Stars: individual: HD 168112 - Stars: 
individual: HD 167971 - X-rays: stars. 


1 INTRODUCTION 

The open cluster NGC6604 is located is the SerOB2 as¬ 
sociation, at a distance of 1.7kpc (iReipurthIl2008fl . It har¬ 
bours two members of the category of particle-accelerating 
colliding-wind binaries (PACWB), i.e. HD 168112 and 
HD 167971. This category of objects includes binary, or 
higher multiplicity systems made of massive stars known to 
be able to accelerate particles up to relativistic energies. The 
mo st up to date catalogue of su ch objects has been published 
by iDe Becker fc Raucd ll2013l l . Even though PACWBs are 
mainly identified through radio observations revealing syn¬ 
chrotron radiation, these objects are potential non-thermal 
emitters in the high energy domain. Relativistic electrons 
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- accele rated via the Diffusive Shock Acceleration (DSA) 
process (lEichler fc Usovl[l9^ : iReitberger et al.l[2014l ) - can 
indeed be responsible for the production of high energy pho¬ 
tons through inverse Compton scattering, with resulting en¬ 
ergies depending on the frequency of t he seed photons and on 
the e nergy of the scattering electrons (IBlumenthal fc Gouldl 
Eiz]). However, it has been demonstrated that the detec¬ 
tion of such a non-thermal soft X-ray emission component is 
very unlikely because of the strength of the th ermal emission 
due to the colliding-winds (iDe Beckeill2007t) . Non-thermal 
X-rays are thus only expected to be detected in the hard 
X-ray do main (above lOkeV ), as demonstr ated in the cases 
of r/ C ar llViotti et al.l l2004l) and WR 140 llSugawara et al.l 

Colliding-wind binaries are indeed well-known ther¬ 
mal X-ray emitters (IStevens et al.l 1 19921 : IPittard fc Parkinl 
I 2 OIOI) . Their X-ray emission comes from the addition of 
the thermal emission from the stellar winds of individual 
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Figure 1. EPIC-MOSl (left) and EPIC-pn (right) images of NGC6604 between 0.3 and 10.0 keV. The positions of the two targets 
are shown, along with the spatial filters used for the extraction of source and background spectra. HD 168112 is located in the upper 
left part, and HD 167971 is the brightest source in the lower right part. Only the September 2014 images are presented. We refer to 
IPe Becker et akl ll2004l . [ioORh for images of the previous epochs. North is up and the East is on the left. 


Table 1. Journal of observations of NGC6604. 


Observ. Id. 

Date 

ID 

Performed duration (s) 

Eff. exp. 

time (s) 

Soft proton flare 




MOS 

pn 

MOS 

pn 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

0008820301 

7 April 2002 

A02 

13120 

10550 

9560 

6820 

yes 

0008820601 

9 September 2002 

S02 

13670 

12050 

13400 

10700 

no 

0740990101 

9 September 2014 

S14 

25500 

23840 

24300 

19000 

yes 


stars with that of the colliding-wind region (CWR). Indi¬ 
vidual winds emit X-rays because of intrinsic shocks re¬ 
sulting from the line-driving instability, leading to slabs 
of win d material characteri zed by temperatures of a few 
10® K llFeldmeier et al.lll997l) . On the other hand, the strong 
shocks due to colliding winds, with pre-shock velocities go¬ 
ing up to 2000-3000 km s“^, are able to ri se post-shock 
plasma temperatures up to several 10^ K llStevens et al.l 
I 1992 I : IPittard Parkinll201Cll ). Both systems investigated 
in this paper, HD 168112 and HD 167971, have been stud- 
ied in X-rays using XMM-N ewton more than a decade ago 
llPe Becker et al]l2004l . [2005l ). Considering their multiple na¬ 
ture with rather long orbital time-scales, a longer term in¬ 
vestigation was thought relevant. For this reason, they were 
observed again about 12 years later to investigate long-term 
trends of their X-ray properties. 

The objective of this study is to complete the descrip¬ 
tion of the X-ray behaviour of HD 168112 and HD 167971 
with respect to previous studies, in ligh t of new results 
published recently about their multiplicity (iDe Becker et al] 

l2012l : [Tbanoglu et aI.l20l3 : ISana et al.l[201^ . As well-studied 

PACWBs, they constitute potential targets for future stud¬ 
ies aiming at modelling the particle acceleration and non- 
thermal emission processes at work in colliding-wind bina¬ 
ries. In this context, X-ray observations constitute a relevant 
probe of the properties of the shocks both responsible for 
their thermal X-ray emission and particle acceleration pro¬ 
cesses. The observations and data processing are described 
in Section[21 The results obtained for both targets are pre¬ 
sented and discussed in Sections[3] and ID The summary and 
conclusions are finally given in Section^ 


2 OBSERVATIONS AND DATA PROCESSING 


NGC6604 has been a target for XMM-Newton at three dif¬ 
ferent dates, two in 2002 and one in 2014. The journal of 
observations is presented Table[T] The observation identifier 
(col. 1) and date (col. 2) are provided, along with an ab¬ 
breviated identifier (col 3) that will be used throughout the 
paper to distinguish the three epochs. Observations A02 and 
S14 were affected by soft proton flares (col. 8), justifying the 
significant difference between the performed duration (col. 
4 and 5) and the effective exposure ti mes (col. 6 and 7) 
of the observations . EPIC instruments (iTurner et aLlliooH : 
IStrfider et al.lf200ll ) were operated in Full Frame mode, and 
the medium filter was used. Data were processed using the 
XMM-Newton Science Analysis Software (SAS) v.14.0.0 on 
the basis of the Observation Data Files (ODF) provided by 
the European Space Agency (ESA). Event lists were filtered 
using standard screening criteria (pattern 12 for MOS and 
pattern ^ 4 for pn). The data processing was performed us¬ 
ing version 3.232 of the Calibration Database (release of 9 
December 2014). Concerning the exposures affected by soft 
proton flares, light curves based on events with Pulse In- 
variant energy cha nnels higher than 10000 were extracted 
dLumb et akl 200^ ). For the April 2002 observation, time 
intervals with background level higher than 25 counts for 
MOS, and 150 counts for pn, were rejected. In September 
2014, the rejection thresholds were 20 counts for MOS and 
150 counts for pn data sets. We note that for the S14 obser¬ 
vation, the aim point was set on the position of HD 167971 
but the signal in RGS spectra was too low to provide data 
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of sufficient quality to warrant adequate spectral analysis. 
This study will therefore focus on EPIC data. 

Our two targets appear in the field of view of the 
three observations, on different locations on the detector. 
In Figure[T] HD 168112 is located in the upper left part, and 
HD 167971 appears in the lower right part. The events re¬ 
lated to the source and background were extracted from the 
complete event list, without applying any energy band filter 
at this stage. 

The source extraction regions for both targets were cir¬ 
cles with 60 arcsec radius, and background spectra were ex¬ 
tracted in annular regions centered on the position of the 
source (with inner radius of 60 arcsec, and outer radius of 
85 arcsec to obtain a background area very similar to that 
of the source). In the case HD 167971, a 12 arcsec radius 
circular region was used to reject a point source located at 
coordinates [18:18:03.089,-12:14:38.90]. For HD 168112, a 30 
arcsec circular region was used to reject a point source at 
coordinates [18:18:44.414,-12:07:53.31] from the background 
region. These regions are shown in Fig.[T]for the EPIC-MOSl 
and EPIC-pn data set, in September 2014. It must however 
be noted that for the EPIC-pn data in September 2014, dif¬ 
ferent background regions were checked. Considering the po¬ 
sition of the source close to the field of view boundary, the 
spectrum is contaminated by a fluorescence line at about 
8keV, part of the so-called internal quiescent background 
(see the XMM-Newton User’s Handbook and iLumb et al.l 
l2002|j . Beside the ’classical’ annular region described above, 
different circular regions were tested to try to get rid of the 
fluorescent feature, but the instrumental feature was still 
present. It is indeed very difficult to make those fluorescent 
features disappear. It must however be noted that the pres¬ 
ence of that feature - albeit visible in the spectrum (see lower 
panel of Fig. [2]- is not significant in the context of the spec¬ 
tral analysis described in Section[3]2] The spectral modelling 
is especially influenced by the real spectral features observed 
at lower energies, and contributing much more significantly 
to the EPIC spectrum. 

The count rates obtained in the three EPIC data sets 
for both targets are quoted in Tabled Please note that the 
count rates do not translate the actual physical fluxes of the 
sources. These numbers are substantially affected by the po¬ 
sition on the detector due to an non-homogeneous response. 
In particular, sources located at the center are more exposed 
than sources closer to the boundary of the field. This has 
to be taken into account while discussing the variability of 
sources located at different detector positions depending on 
epoch. For this reason, the discussion about the long term 
behaviour (over about a decade) of the targets in X-rays 
will be based on the results of the spectral modelling (see 
Sections l3.3.2l and 14. 3. 2D . In addition, one has to be espe¬ 
cially careful for sources located close to the boundaries of 
the field of view. For instance, it is clear in Figure[T] that 
HD 168112 is not optimally exposed in September 2014 for 
MOS instruments. This leads to a rather weak signal in the 
harder part of the EPIC bandpass where the source is less 
bright, as shown in spectra presented in the next sections. 

For all spectra at every epoch, response matrices and 
ancillary response files were computed using the rmfgen 
and arfgen metatasks. Spectra were grouped to get at least 
20 events per energy bin. The spectral ana lysis was per- 
formed using the XSPEC software (v.12.8.2. lArnaudlll996l : 


Table 2. Count rates measured by the three EPIC instruments 
(between 0.3 and lOkeV) for the two targets at the three obser¬ 
vation epochs. 



HD 168112 

HD 167971 

April 2002 

MOSl 

MOS2 

pn 

0.059±0.003 cnts-i 
0.055±0.003 cnts-i 
0.173±0.007 cnts-i 

0.148±0.004 cnts-l 
0.144±0.004 cnts-l 
0.401±0.010 cnts-l 

September 2002 

MOSl 

MOS2 

pn 

0.036±0.002 cnts-i 
0.041±0.002 cnts-l 
0.126±0.005 cnts-l 

0.138±0.004 cnts-l 
0.160±0.004 cnts-l 
0.476±0.008 cnts-i 

September 2014 

MOSl 

MOS2 

pn 

0.034±0.002 cnts-l 
0.035±0.002 cnts-i 
0.106±0.003 cnts-i 

0.245±0.004 cnts-l 
0.235±0.003 cnts-i 
0.834±0.008 cnts-i 


iDorman fc Arnaudll200l] f . Spectral fitting was executed us¬ 
ing the Xv value as goodness-of-fit indicator, which is ade¬ 
quate as the grouping criterion leads to a sufficient number 
of counts in each energy bin to assume a Gaussian distri¬ 
bution. All error bars were computed on the basis of the 
error command within XSPEC. This leads to an adequate 
estimate of confidence intervals of model parameters based 
on a systematic exploration of the parameter space. Such a 
procedure allows to refine the spectral fitting and the con¬ 
vergence to a best-fit model. Details on the spectral analysis 
for HD 168112 and HD 167971 are provided in Sections 13.21 
and 14.21 respectively. 

3 HD 168112 
3.1 The system 

HD 168112 (BD-12 4988) includes an 05.5111(1+) primary 
with an unidentified secondary. This star is known as a non- 
thermal radio emitter since the large campaigns in the 1980’s 
dedicated to massive stars (iBieeing et al.lll989l ). with a well- 
established negative spectral index pointing to a clear syn¬ 
chrotron signature. In the context of the so-called ’standard 
model’ for particle acceleration, this object was suspected 
to be a binary with a non-thermal emission region coinci¬ 
dent with a colliding-wind region. However, t he search for 
a spe ctroscopic companion failed to detect any llRauw et al.l 
I 2 OO 5 II . The solution came from another technique with the 
development of long baseline interferometry which allowed 
to resolve a companion in HD 168112 usi ng the Very Large 
Telescope Interferometer (VLTI, Chile). ISana et al.l (l2014l l 
reported on the presence of a companion with a separation 
of a 3.33±0.17 mas in May 2012, providing evidence for the 
binary nature of this object. At a distance of 1.7 kpc, this 
separation translates into a projected linear separation of 
the order of 6 au. The interferometric study revealed a mag¬ 
nitude difference AH = 0.17 ± 0.19. On the basis of the syn- 
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thetic photometry published bv iMartins fc Pl^ (l2006l i. this 
confidence interval of AH suggests spectral classifications for 
the secondary ranging between 05.5III and 07.5III. One ex¬ 
treme value of AH is also compatible with an 04V spectral 
type for the secondary, but such a classification would have 
been revealed by the existing visible spectra of the system, 
but this is not the case. 

Before the advent of VLTI results, additional hints for 
the binary nature came from XMM-Newton observation re¬ 
vealing a spectrum too hard to be compatible with that 
of a single O-type star, along with a weak variation of 
the X-ray flux between two epochs separated by 5 months 
toe Becker et al]l2004^ ■ Such an intrinsic variation from an 
individual stellar wind is indeed very unlikely. On the other 
hand, radio observat ions revealed also a significant vari- 
ability of HD 168112 toe Becker et al]l2004l : iBlomme et al.l 
l2005h . On the basis of the r adio light curve explo iting data 
spread over about 20 years, iBlomme et al.l ll2005h explored 
potential orbital time-scales but the time series was too 
sparse to be conclusive. So far, the orbital period is still 
undetermined. 

When one is dealing with colliding-winds in massive bi¬ 
naries, i t is relevant to introd uce the cooling parameter de¬ 
fined bv IStevens et al.l (Il992l i. This parameter is basically a 
measurement of the ratio of the cooling time of the shocked 
gas over the flow time. A value lower or close to one points 
to a radiative regime where the X-ray emission is typically 
a fraction of the kinetic power injected in the wind collision 
(whatever the stellar separation). Alternatively, a cooling 
parameter significantly larger than one characterizes an adi¬ 
abatic regime, with an X-ray emission from the shocked gas 
inversely proportional to the stellar separation (the so-called 
1/D-law, where D is the separation). Considering the typi¬ 
cal values for the mass loss rate and terminal velocity of the 
05.5HI primary, on e gets M = 1.35 x 10~^Mf7iyr~^, and 
Voo = 2800 km s“^ dMuiires et al. 2012ll. Th e definition of 
the cooling parameter given by Stevens et al. 1(H) can be 
translated into the following criterion that gives a lower limit 
on the distance (in cm) between a star and the stagnation 
point warranting an adiabatic regime: 


10^1 


M 


V 


4 

oo 


with quantities expressed in the same units as above. This 
translates into a minimum distance of about 0.015 au. Even 
considering that in a highly eccentric orbit the winds may 
collide at a very short distance from the star with the weaker 
wind at periastron, dmin is much smaller than the projected 
separation measured with the VLTI, suggesting an adia¬ 
batic regime. One is lacking the required accurate informa¬ 
tion to apply this criterion to the secondary object. How¬ 
ever, assuming the spectral type range mentioned above 
in this Section, the wind parameters should be similar to 
that of the the primary or take the following values if the 
spectral type is 07.5HI: M = 2.5 x 10~^ Mq yr“^ and Voo 
= 4400kms“^ llMuiires et al.l[2012l l. The latter assumption 
leads to a lower limit on the distance between the colliding- 
wind region and the secondary of about 0.0004 au for the 
adiabatic regime (much smaller th an the projected separa¬ 
tion measured bv ISana et al.|[^14l '). One may therefore rea¬ 
sonably consider that the expected properties of the shocked 


winds in HD 168112 point to the adiabatic regime even in a 
highly eccentric system. 


3.2 X-ray properties 


It is worth to point to some differences with respect to the 
analysis performed bv iDe Becker et al.l (120041') on the basis 
of the A02 and S02 data sets. Notably, the event lists were 
filtered to get rid of time intervals contaminated by a soft 
proton flare. This leads the A02 exposure to be reduced by 
about 30% with respect to previous analysis (see TableO, 
but the spectra are expected to be cleaner. 

Spectra were first fitted using two-temperature models 
of the type wabs*wind* (apec+apec). The wabs component 
stands for the photoelectric absorption of X-rays by the in¬ 
terstellar mater ial along the line of sight, and it is based on 
abundances bv lAnders fc Eb ihara (Il982|). U sing t he same 
approach as described by De Becker et al.l (|2004| L it was 
fixed to a value of 0.58 x 10^^ cm“^. This value was obtained 
assuming Nh = 5.8 x 10^^ x E(B - V)cm“^ (iBohlin et al.l 
, with a col our excess determined f rom an observed 
(B - V) of -1-0.69 (IChlebow ski et a l]ll989l~) and an intrinsic 
colour (B - V)o of -0.309 ( Mihalas fc Bimi^ 1 198 il l. The 
wind component accounts for the absorption of X-rays by 
the wind (i.e. circumstellar) material. This absorption com¬ 
ponent was computed assuming the presence of ionized ma¬ 
terial made of the ten most abundan t elements with solar 
abundances (lAnders fc Grevesselll98^. Details on this ion - 
ized absorption model are given in De Becker et al.l (l2005ll . 
and references therein. The apec emission components are 
used to reproduce the thermal X-ray emission spectrum from 
an optically thin plasma heated by shocks, either intrinsic 
or due to colliding-winds. A typical plasma temperature of 
about 3 X 10® K is obtained for the soft emission component. 
The temperature of the hard component is less well deter¬ 
mined, with values in the rang e of 15-30 x 10® K. We note 
that in the approach followed bv iDe Becker et"^ (l2004l i the 
fit were not optimized by an exploration of the parameter 
space using the error command within XSPEC. The explo¬ 
ration of the best-fit parameter space revealed convergence 
difficulties which were not apparent at the time of the pre¬ 
vious analysis. The problem comes from the existence of dif¬ 
ferent local minima, in particular characterized by kT values 
of about 0.3 or 0.7 keV for the soft component. This happens 
only for a few spectra in the series. Typically, for the lower 
quality spectra, a degeneracy may appear between models 
with different temperatures weighted by different normal¬ 
ization parameters. The kT value of 0.3 keV appears to be 
consistently adequate for all spectra at every epoch. How¬ 
ever, the highest quality spectra did not stabilize in solutions 
including a plasma temperature of about 0.7 keV for the soft 
component. For these reasons, the solution with the lower 
value for k Ti was privileged, in agr eement with the results 
obtained bv iDe Becker et al.l ((200J). 

Best-fit parameters are summarized in Table[3] for 
EPIC-MOSl and EPIC-MOS2 data simultaneously, for 
EPIC-pn individually and for the combined EPIC instru¬ 
ments. Fluxes have also been determined along with error 
bars based on the exploration of the parameter confidence 
ranges of the best fit. Globally, the soft plasma temperature 
is constant throughout the time series. The hard tempera¬ 
ture does not seem to present a clear trend as a function of 
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Table 3. Best-fit parameters for HD 168112 using the wabs*wind*(apec+apec) model (in the 0.3-10 keV energy band). The error bars 
represent the 90 % confidence level. 


Instr. 

(10^2 cm-2) 

kTi 

(keV) 

Ni 

(10-®cm-5) 

kTa 

(keV) 

N2 

(10-®cm-5) 

xl (d.o.f.) 

fobs 

(ergcm”^ 

April 2002 

MOS 

o.7to°i® 

0.291006 

3 i50+^ ’^O 
o.ou _2 67 

3.68111® 

0 31+0.10 
^•'J-L_o.09 

1.10 (59) 

5.4l®'3 x 10-1® 

pn 

0.81+01^ 

0.2810-06 

8.13l®+® 

2.471J;?® 

0 53+0-14 
U.OO_Q 17 

1.08 (70) 

6.51°;® X 10-1° 

EPIC 


0.2910-03 

6.551+1° 

3 34+2.12 

0 34+0.08 

1.15 (142) 

5.81°;^ X 10-1° 

September 2002 

MOS 

0.53t°;l^ 

0 30+0 04 
U.OU_Q Q 2 

2 60"^^'^*^ 
^•ou_Q 35 

2 541-94 
^•O^_0.77 

0 26^*^'*^^ 
0 '-^O_Q Q 7 

0.97 (55) 

3.8l°;g X 10-1° 

pn 

o.8to°®i 

_0.04 

7 rq+12.00 

1 •oo_i 12 

z.yo_Q 93 

n 00+0.12 

0 -oo_Q Q 9 

0.77 (76) 

5.612;° 10-1° 

EPIC 

0.72+012 

^•^O_o.o2 

‘ -^^-3.82 

2 45+0-93 

0 33+0.07 

1.07 (149) 

4.4l° ° X 10-1° 

September 2014 

MOS 

'-'•^^-0.18 

fi qo + 0.03 

u.oz_0 03 

9 i::q + 1-26 
z.oo_0 92 

9 11 +0.52 
^•^-^-0.25 

0 

l^-'^-O.lO 

1.02 (85) 

7.3l°;^ X 10-1° 

pn 

^•^^-0.16 


2.6411®® 

1 ^7+0.14 

-0.13 

n 7 C+O .07 

'O_o.o8 

1.15 (115) 

5.91°;® X 10-1° 

EPIC 

14 

0 31+0.04 

9 irq+1.39 

i.yu_o 17 

r\ 79 + 0.06 

1^- ''^_0.05 

1.15 (218) 

6.6l°;i X 10-1° 


time, as discrepancies appear between different data sets of a 
given epoch. The most plausible variations concern the wind 
absorption component and the normalization parameters 
of the emission components. On average, the faintest and 
brightest spectra are associated to the S02 and S14 epochs, 
respectively. EPIC spectra with the associated best-fit mod¬ 
els are shown in Fig. [2] The use of a three-temperature model 
did not improve the results of the spectral fitting. Typi¬ 
cal reduced values were indeed slightly lower than those 
quoted in Tabled For instance, for combined EPIC instru¬ 
ments the values were 1.15, 0.98 and 1.03 respectively for 
the three epochs. However, the three emission components 
model led to kT values of about 0.2-0.3 and 0.8-1.0keV, 
with a very faint third component characterized by an ex¬ 
cessive kT value (8 - 15 keV) with huge relative errors (often 
close to 100%). This suggests clearly that the use of a third 
component does not provide any improvement of the physi¬ 
cal inter 2 retetioiirf the present data set. Finally, as justified 
in IPe Becker! ll2007l) . the use of a power law component to 
account for a potential inverse Compton scattering process 
was not considered as one is dealing only with soft X-rays 
(below 10 keV) dominated by thermal emission. 


3.3 Discussion 

3.3.1 The X-ray luminosity 


1.7 X 10®®ergs“^. As the spectral classification of the sec¬ 
ondary is not determined, one can only speculate on its 
main properties. On the basis of the discussion in Sec- 
tion l3.ll a likely range for the secondary bolometric lumi¬ 
nosity is 0.9-1.7 X lO^^ergs”^. These values translate into 
a Lx/Lboi ratio of about 2.6-3.8 x 10~^, pointing a to a 
significant - although weak - X-ray overluminosity, with 
respect to the L x/Lbof ratio of ^ 10~^ expected for sin¬ 
gle O - star winds dChlebowskill 19891 : IChlebowski fc Garmanvl 


[l99ll : lOwocki et al. 2013h . This results demonstrates that 
the colliding-wind region in the HD 168112 system is mod¬ 
erately bright in X-rays. This is in agreement with the global 
shape of the EPIC spectrum which presents a noticeable - 
albeit not spectacular - hardness. 

The soft X-ray emission (between 0.3 and 10.0 keV) 
from the colliding winds could be roughly estimated, sub¬ 
tracting the expected contributions from individual stars. 
Assuming a 10“^ luminosity ratio for the stars in the sys¬ 
tem, we obtain 1.7 x lO^^ergs”*^ for the primary and 0.9- 
1.7 X 10^^ ergs“^ for the secondary. These numbers lead to 
an X-ray luminosity for the colliding-wind region of about 
6.5 X 10®^ ergs“’^. 

The latter number could be compared to the kinetic 
power of the stellar winds: 


The overall X-ray emission from HD 168112 corrected 
for interstellar absorption corresponds to a flux of 
2.6-3.0 X 10“^^ ergcm“^ s”*^. Considering the distance to 
NGC 6604 is 1700 pc, this translates into a luminos¬ 
ity (Lx) range of 0.9-1.0 x 10®® ergs“^. On the basis 
of the calibration of parameters given by iMartins et al.l 
ll2005l) . the bolometric luminosity of the primary should be 


This quantity can be calcul ated using the typic al stellar 
wind parameters published bv iMuiires et al.l (l2012l) . as used 
in Section ISTi One obtains a total (for both stars) kinetic 
power in the range of 4.8-6.6 x 10®® ergs“^, translating into 
a Lx/Pfein ratio of 0.9-1.5 x 10““*. These numbers are in 
agre ement with the expected o rder of magnitude mentioned 
by toe Becker fc Raucd l2013l ). Such a ratio can not be 
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Energy (keV) 



Energy (keV) 

Figure 2. Fit of the two-temperature model to EPIC-pn spectra 
obtained in September 2002 (top) and September 2014 (bottom) 
between 0.3 and 10.0 keV. In both cases, the lower panel shows 
the residuals in the sense data minus model. The apparent spec¬ 
tral feature at about 8 keV in September 2014 results from the 
fluorescence of the pn enclosure hit by charged particles. 


used as a immediate proxy for the power injection into the 
colliding-winds, as it will depend notably on the geometry of 
colliding-wind region which is not determined in this poorly 
characterized system. However, the expected future deter¬ 
mination of the orbital parameters of HD 168112 along with 
ongoing theoretical developments should allow us to address 
with more details the issue of the energy budget in this sys¬ 
tem. The latter point is crucial the discuss the energy in¬ 
jection rate into non-thermal processes, and in particular in 
the particle acceleration process which takes place in this 
system. 

3.3.2 Variability 

Beside the main spectral properties (essentially in terms of 
plasma temp erature) found in agree ment with the results 
published bv iDe Becker et al.l (l2004l) . the results presented 
in Section[3]2] deserve to be discussed from the variability 
point of view. In the case of a binary system, the most 
straightforward source of variability could be attributed 
to the orbital motion. In an eccentric orbit, the phase- 
dependent separation will lead to changes in the physical 
properties of the colliding-wind region where a significant 
fraction of the X-ray emission is produced. As explained in 


Section ISTi in the adiabatic regime the intrinsic X-ray emis¬ 
sion (corrected for both interstellar and wind absorption) 
is expected to be inversely proportional to the stellar sep¬ 
aration. On the other hand, the absorbing column made of 
wind material is expected to change with the orbital phase. 
This orientation effect will also induce variations in the ob¬ 
served spectrum. However, provided that pre-shock veloc¬ 
ities of the colliding winds do not change, one should not 
expect a change in the plasma temperatures from one epoch 
to the other. In an eccentric orbit, this must be true provided 
the winds have reached their terminal velocities before col¬ 
liding at all observed orbital phases. 

As a consequence, the clearer trends are provided by 
the wind absorption component and the normalization (re¬ 
lated to the intrinsic X-ray luminosity) of the emission com- 
ponent(s). The former point is clearly an orientation effect 
which depends intimately on the geometry of the system as 
seen by the observer. Without any indication on the orbital 
parameters of HD 168112, the interpretation of these trends 
in terms of geometrical consideration would be highly spec¬ 
ulative. 

Concerning the varying intrinsic X-ray luminosity, the 
values can be estimated on the basis the best-fit parameters 
for EPIC instruments found in Table|3] Intrinsic X-ray fluxes 
are 1.49 x 10“'^^, 1.51 x 10“^^ and 0.68 x 10“^^ ergcm“^ s“^ 
for the A02, S02 and S14 epochs. For a distance of about 
1700 pc, these fluxes translate into X-ray luminosities of 
5.2 X 10^®, 5.2 X 10^® and 2.4 x 10^® erg s“’^. In the context of 
adiabatic shocked winds this quantity should follow a 1/D 
trend. It is therefore reassuring to see that the two 2002 
observations give an almost identical value of the intrinsic 
X-ray luminosity, and the significantly different value for 
2014 points very likely to a significantly different separation 
between the two stars at these two epochs. 

The eccentricity (e) of the orbit is related to the sep¬ 
aration ratio at apastron (Da) and periastron (Dp) by the 
following identify, 

/ Dp \ _ 1 — e 

\Da/ r 1 T e 

where the index r stands for ’real’, as opposed to m for 
’measured’ (see below). 

On the other hand, on the basis of the 1/D-law, one can 
write that the measured separation ratio is 



where Lx,i4 and Lx,02 are the minimum and maximum mea¬ 
sured intrinsic X-ray luminosities. 

Considering that 


f —1 

VDaA 




we have 


e ^ 


^-1 

^+1 


With the ^ value estimated above, one obtains a rough lower 
limit on the eccentricity equal to 0.38. 

One should also mention that the measured X-ray emis¬ 
sion does not come only from the colliding-wind region, but 
individual winds contribute also. The variable part of the 
overall emission appears thus on top of a constant emission. 
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It means that the ^ value constitutes itself a lower limit on 
the flux ratio coming from the colliding-wind region alone. 
The above lower limit on the eccentricity is therefore highly 
conservative. 


4 HD 167971 
4.1 The system 

HD 167971 (BD -12 4980, MY Ser) has been known as a 
po tential triple sy s tem s ince the 1980’s, with t h e stu dies 
by iLeitherer et al.l (Il987l l and iDavidee fc Forbe^ (Il988l ~) . It 
consists in a close eclipsing binary with a period of about 
3.32 days, with a third object evolving on a much wider 
orbit. First evidence for the gravitational link between the 
close pair and the third star was provided by VLTI obser¬ 
vations revealing a t ime-dependent astrom etry suggestive of 
an orbital mot ion Becker et al.l l2012l l. More recently, 
llbanoglu et al.l ll2013l ) published the very first orbital solu¬ 
tion for the wide orbit. According to this spectroscopic solu¬ 
tion the period should be 21.2±0.7yr, with an eccentricity 
of 0.53±0.05 and a projected semi-major axis of 5.5±0.7 au. 
The spectral types are 07.5III and 09.5III for the close 
pair, and the spectral classification of the third component 
is O9.5-B0(I or III). _ _ 

The orbital period determined bv llbanoglu et aP ll2013l l 
is in excellent agreement with the variation time-scale sug¬ 
gested by the radio light curve published bv iBlomme et aP 
(1200711 . This lends strong support to the idea that the non- 
thermal radio emission occurring in the wide orbit. In such a 
hierarchize d system, one may disti nguish two colliding-wind 
regions fsee lPe Becker et al.ll2013 for a discussion). The sit¬ 
uation is illustrated in Fig.|3] CWab and CWAB stand for 
the two expected colliding-wind regions, in the short and 
long period orbits respectively. 

The analysis of the X-ray spectra a t epochs A02 and 
S02 performed bv iDe Becker et al.l ll2005l l revealed a signifi¬ 
cantly hard spectrum, with a warm component temperature 
of several 10^ K, in addition to a dominant softer emission. 
It suggested also a slight - but marginal - flux variation 
between the two epochs. The orig in of that variat i on re - 
mained unclear. At the time of the lPe Becker et al.l (l2005l l 
study, no orbital parameter was known for the third object. 
It was not even clear that there was a gravitational relation 
between components A and B. Even though a possible sce¬ 
nario could be an orbital modulation of the X-ray emission 
from CWAB (see Figure[3|, a potential partial eclipse within 
the close Aab pair could not completely be rejected. 

As far as colliding-wind regions are involved in the X- 
ray emission in HD 167971, the issue of the regime of the 
shocked wind - radiative or adiabatic - is worth a discussion. 
Within the close Aab pair, the cooling parameter defined by 
IStevens et al. I (Il992ll can easily be estimated. The typical 
wind parameters for c ompo nents Aa and Ab can be taken 
from iMuiires et al.l (l2012l l: MAa = 2.5 x 10“^MQyr“^, 
Mas = 1.4 x 10“^MQyr“^, Voo,Aa = 4400kms“^ and 
Voo.Ab = 2600 kms“^. According to llbanoglu et al.l ll2013ll . 
the projected semi-major axis of the eclipsing binary (with 
a circular orbit) is about 36 R0. Such a short distance does 
certainly not allow the winds to reach their terminal veloci¬ 
ties before colliding. The pre-shock velocities must therefore 



Figure 3. Schematic view of the triple system HD 167971. The 
thick black curves represent the shocks of the collidlng-wind re¬ 
gions. Two colliding-wind regions are present: CWab between 
components Aa and Ab, and CWAB between components A and 
B. The relative stellar separations and sizes are not to scale. 


be estimated on the basis of the classic al /1-law for win d 
velocities, assuming /3 = 0.8 for 0-stars (|Puls et al.l l2008l 'l: 
i.e. VAa = 1300kms“^ and Vas = 900kms“^. These values 
lead to cooling parameter values of the order of unity, point¬ 
ing to a radiative regime. Concerning CWAB involving the 
shocked wind of component B, the much larger separation 
leads inevitably to an adiabatic regime. 


4.2 X-ray properties 


EPIC spectra were analyzed according to the same approach 
as described in Section[33] Good results were obtained us¬ 
ing a two-temperature model made of the same components 
as for HD 168112. For the insterstellar absorption, the same 
approach as described in Sect. l3.2l was followed. On the ba - 
sis an observed (B - V) of -1-0.77 (IChlebowski et ah 1989|) 
and a n intrinsic colour (B - V)o of-0.31 ( Mihalas fc Binn^ 
Il98lll . a Njy value of 0.63 x 10^^ cm“^ was determined and 
adopted for all models (see lPe Becker et al.ll2005l l. The in¬ 
terstellar column was fixed to that value, allowing all other 
parameters to vary. The results are summarized in TableH] 

The best-fit temperatures (kT) are about 0.3 keV and 
1.0-1.2keV, respectively for the soft and hard components. 
The consistency between results obtained with MOS and pn 
data is very good, except for the S02 epoch. The measured 
flux is very similar at all epochs, with only a slight decrease 
in S02. However, this trend is present only in MOS data, 
and not in the pn data set. It must be reminded that in the 
IPe Becker et al.l (l2005l l study, only MOS data were consid¬ 
ered for epoch A02 and only pn data were considered for 
epoch S02. The reason could be the poor handling of gaps 
and bad columns be the first releases of the Science Analy¬ 
sis Software. Even though the data processing with the most 
recent releases of the SAS is expected to apply an adequate 
correction, the quality of the flux measurement may be ques¬ 
tioned. For this reason, we will refrain to base our discussion 
on the pn result for the S02 epoch. 

The use of a three-temperature model did not improve 
the results on the complete time series. The typical plasma 
temperature of the soft and medium components were con¬ 
sistently of the order of 0.25-0.30 keV and 0.95-1.00 keV, 
respectively. However, the lack of significant signal in the 
higher energy part of the EPIC spectrum led to a difficulty 
to constrain both the temperature and the normalization 
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Table 4. Best-fit parameters for HD 167971 using the wabs*wind*(apec+apec) model (in the 0.3-10 keV energy band). The error bars 
represent the 90 % confidence level. 


Instr. 

"^wind 
(10^2 cm-2) 

kTi 

(keV) 

Ni 

(10”® cm“" 

kTz 

') (keV) 

N2 

(10”® cm-®) 

xi (d.o.f.) 

fobs 
(ergcm“^ 

April 2002 

MOS 

u-0o_o 10 


-| qq+1.07 

±.oo_o 54 

0 99"’'*^-^*^ 
'^•^^-0.07 

2 43+0-41 

1.14 (121) 

1 - 6 i 1®;?6 ^ 10 "^^ 

pn 


n 9'7+0.03 

-0.03 

1 04+0.77 

1 16+0-10 
-*■•^^-0.22 

9 97+0.78 
-0.28 

1.09 (129) 

1.681°;°^ x 10-1® 

EPIC 

0 5R+O.07 

02 

1 41+0-66 

1 00'*"*^'^^ 

9 1=7+0.27 
-0.31 

1.16 (266) 

1.62l°;0® x 10-1® 

September 2002 

MOS 


U.DU_o 03 

0 92'^*^'^® 
'^•^^-0.30 

1.201°®® 

1 74+0.20 
!• '^-0.19 

1.10 (162) 

l.43l®;°0 ^ 10"^^ 

pn 

0 47+0 08 
0'-^<_0.o8 

0.25t°®® 

1 -17+0.44 
-*-•-*-'-0.34 

0-98l°®^ 

i-yo-0.13 

0.93 (220) 

1.34l®;°® x 10-1® 

EPIC 

0 49+0 06 

0.28+°®® 

-^•^^-0.23 

-1 1 9+0.09 

-*-•-*-^-0.15 

'^-0.13 

1.09 (395) 

1.39l°;°® x 10-1® 

September 2014 

MOS 

0'-O^_0.05 

n 97'+0.02 
'^•^'-0.01 

1 qc:+0.28 

-L-'9d_o 23 

-10+0.04 

l-lo_o 04 

2-ool°:®i; 

1.41 (259) 

1.57lH^ X 10-^^ 

pn 

0 59+0 04 

0.28+°®® 

1 q9+0.24 

-^•'^^-0.23 

-1 99+0.04 

-*-•^^-0.04 

2.05l°®® 

1.08 (402) 

l.6ll°;°j x 10-1® 

EPIC 

U.O-i_o.o4 

0-28l°®®i 

1 97+0.17 
-*-•^'-0.17 

i-^oi_o.o3 

2-oil2:S^ 

1.38 (685) 

l.59l°;°® X 10-1® 


parameter of the hardest component. The only relevant re¬ 
sult with a three-temperature model was certainly obtained 
for the September 2014 data set, where HD 167971 is well 
exposed at the center of the field of view. For this specific 
data set, some improvement was obtained by switching to a 
three-temperature model. The results are presented in Ta- 
ble[5l It is notably noticeable in the plot shown in Figure|4] 
that the iron K line at about 6.7 keV is well identified for this 
data set. That feature in not so obvious in other spectra of 
the time series, certainly explaining the difficulty to fit the 
hardest thermal component. For the same reason as given 
in the case of HD 168112, no non-thermal emission compo¬ 
nent was used to investigate any potential X-ray emission of 
inverse Compton scattering origin. 

EPIC spectra obtained at the three epochs are shown in 
Fig-El The wide orbit of t he triple system calcula ted on the 
basis of the ephemeris of llbanoglu et aP (l2013l) is shown, 
along with the corresponding orbital phases of the obser¬ 
vations (see Sect. 14.3^ for a discussion). The dashed lines 
represent the stellar separation between components A an 
B, emphasizing the significant changes between 2002 and 
2014. 


4.3 Discussion 

4-3.1 The X-ray luminosity 

Consi dering the spectral types given by llbanoglu et al.l 
1I2OI3I) for the components of HD 167971, typical bolometric 
lu minosities can be estim ated on the basis of the calibration 
bv iMartins et al] (I2OO5I ). These values, collected in Table|6l 
lead to a bolometric luminosity of 1.9 - 2.6 x 10®® erg s“^. On 
the other hand, the X-ray flux (in the 0.3-10.0 keV range) 
corrected for interstellar absorption can be estimated on the 


Table 5. Best-fit parameters for the September 2014 EPIC- 
pn data set of HD 167971 using a three-temperature model 
(wabs*wind*(ap 6 c-^apec•^apec), in the 0.3—lOkeV energy band). 
The error bars represent the 90 % confidence level. 


Parameter 

Value 


0-45+o°g X 10®® cm-® 

kTi 

0.29t°°®keV 

Ni 

8.26+j ®g X 10-®cm-® 

kT 2 

0.94tO°®keV 

N 2 

1.67+° ®® X 10-®cm-5 

kTs 

1.98l°°®keV 

N 3 

0 - 6 + 0.24 X 10-®cm-5 

xi (d.o.f.) 
fobs 

0.95 (400) 

l-68j:°°® X 10-1® erg cm-® 3-1 


basis of the best-fit parameters, i.e. lO”®*^ ergcm“® s“®. 
For a distance to NGC6604 of 1700 pc, one obtains an X- 
ray luminosity of about 3.5 x 10®®ergs“^. As a result, one 
obtains a Tx/Thai ratio of 1.3 - 1.8 x 10“®, clearly pointing 
to a significant overluminosity certainly due to the colliding- 
wind regions. 

On the basis of the wind parameters mentioned in Sec- 
tion l4.1l for components Aa and Ab, one can estimate a total 
kinetic power for component Aab of about 1.8 x 10®® ergs”®. 
For component B, assuming Mb = 1.4 — 3 x 10”®" Mq yr“® 
and Voo.B = 2600-2900 km s“® depending on the luminos- 
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Figure 4. Wide orbit of HD 167971 based on the parameters published bv llbanoglu et ^ ll2013h . The plane of the page/screen coincides 
with the orbital plane. Periastron passage corresponds to an orbital phase equal to zero and dotted straight lines represent the stellar 
separation at the three epochs. Symbols (triangle, square and hexagon) represent the position of component B on the orbit at the three 
epochs (component Aab is at the focus on the right). The dark arrows point to the corresponding EPIC (MOS and pn) spectra along 
with their best-fit model. Individual model components are represented by the dotted curves. All models are two-temperature models, 
except for the pn spectra in September 2014 which is significantly best represented by a three-temperature model. 
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Table 6. Adopted parameters for the components in HD 167971. 
The two values for the total bolometric luminosity come from the 
two hypotheses for the luminosity class of component B. 


Comp. 

Sp. type 

^bol 

(ergs-1) 

Aa 

07.5III 

8.9 x 10^® 

Ab 

09.5III 

5.1 X 

B 

09.5III 

5.1 X lO^* 


09.51 

12.1 X lO®* 

Total 


19.1 X 10^® 



26.1 X 10®® 


ity class (III or I), we obtain 3 - 8 x 10®®ergs“^. 

The total kinetic power in the triple system is thus ~ 2.1- 
2.6 X 10®®ergs“^. The expected intrinsic emission from the 
individual stellar winds must be about 10“^ times the bolo¬ 
metric luminosities. On the basis of the values quoted in 
Table[6l the expected intrinsic X-ray luminosity is in the 
range 1.9-2.6 x 10®^ ergs”^. Subtracting this contribution 
from the total measured X-ray luminosity corrected for in¬ 
terstellar absorption, the contribution from the colliding- 
wind region is about 3.2-3.3 x 10^® ergs“^. The Lx/Pfein 
ratio is thus in the range 1.2-1.6 x 10~®. It appears that 
the conversion rate of kinetic power into thermal soft X-ray 
emission from the colliding winds is an order of magnitude 
larger than for HD 168112. This can be explained by the fact 
that two colliding-wind regions are present, and in partic¬ 
ular one of them occurs in a close system in the radiative 
regime, highly efficient at producing X-rays. However, it is 
not possible with the present information to disentangle con¬ 
tributions of the X-ray luminosity coming from CWab and 
CWAB. 

As a comparison, Lx/Pfcin ratio can also be estimated 
in the cases of two 0-1-0 PACWBs, one rather short 
period and one long period system. As a short period 
system, let us consider Cyg OB2 #8A (0 6If + 05.5111(f)) 
whos e orbital period is about 22 days ( De^^ecker_jt_^ 
[2OOJ). The X-ray luminosity measured by IPe Becker et al. 
(1200611 was about 1.0-1.5 x 10^"^ergs“^ (overluminous by a 
factor 20-30). For these spectral type s, the estimated ki - 
netic power (based on parameters by iMuiires et al.l l2012l l 
is about 8.5 x 10®®ergs“^. This leads to a Lx/Pfein ratio 
of about 1.2-1.8 x 10~®. For a long period system, the case 
of HD 167974 might be relevant. It consists of an 03.5V -I- 
05-5 .5 V system, w i th a p eriod of about 8.6 yr (iRauw et al.l 
|2012|) . [ilauw et ^ (120021) measured an X-ray luminosity of 
1.4 x 10^^ ergs“^. The spectral types of the components sug¬ 
gest a bolometric luminosity of about 3.4 x 10^® ergs“^. Re¬ 
moving the expected contributions from individual winds 
(10“^ times hhoi), the X-ray contribution from the collid¬ 
ing winds is about 1.1 x lO^^ergs”^. These number allow 
to determine a Lx/Pkin ratio of about 1.0 x 10“^. Among 
these two examples, HD 167971 is closer to the case of 
Cyg OB2 #8A characterized by a much high conversion rate 
of kinetic power to thermal X-rays produced by the colliding 
winds. This emphasizes further the role played by the CWab 
region in the thermal X-ray spectrum of HD 167971. 


4-3.2 Variability 

On th e basis of the ephemeris published bv llbanoglu et ahl 
(l2013li . one can compute approximate orbital phases of the 
long orbit for the three epochs of XMM-Newton observa¬ 
tions. Considering the large uncertainties on the orbital pa¬ 
rameters, it is more realistic to give ranges of values, rather 
than unique values for the orbital phase. These ranges are 
0.57-0.60, 0.59-0.62 and 0.16-0.18, respectively for A02, S02 
and S14. With an eccentricity of about 0.53, the relative sep¬ 
aration between components A and B is expected to drop 
by a factor 1/3 between 2002 and 2014 (see Fig. |4] for a plot 
of the orbit). As a result, the thermal X-ray emission from 
CWAB is expected to be higher in 2014. However, it must be 
emphasized that the total thermal X-ray emission measured 
in the EPIC spectra comes from 5 contributions: (i) the wind 
of component Aa, (ii) the wind of components Ab, (iii) the 
wind of component B, (iv) the wind interaction CWab, and 
finally (v) the wind interaction CWAB. Among these contri¬ 
butions, only the last one is likely to change as a function of 
the orbital phase of the wide orbit. The anticipated change 
in the total measured X-ray flux should therefore be diluted. 
The measured fluxes do not point to any significant change 
of the X-ray flux between April 2002 and September 2014. 

The latter result appears quite surprising consider¬ 
ing the X-ray brightness of the colliding-wind regions in 
HD 167971 suggested by the Lx/L6oi ratio. As a tentative 
explanation, one should note that an important fraction of 
the X-ray overluminosity may come from CWab, with only 
a moderate/minor contribution from CWAB. The short pe¬ 
riod (about 3.32 d) in the Aab pair favours indeed a signif¬ 
icantly high emission measure in CWab, in addition to the 
fact that the radiative r egime is more efficient at producing 
therm al X-rays (see e.g. IStevens et al.l[l9^ : lAntokhin et al.1 
l2004h . Considering the short separation in the eclipsing bi¬ 
nary, the stellar winds have no chance to have reached their 
terminal velocity (see Section lTTII) . The post-shock temper¬ 
ature is therefore lower than in wider systems. The thermal 
X-ray emission from CWab is therefore expected to be rather 
soft. The fact that the global EPIC spectrum of HD 167971 
is not so hard supports the idea that the harder component 
expected to come from CWAB (with higher pre-shock ve¬ 
locities) does not dominate the X-ray emission. The plasma 
temperatures reported in Table|4]are indeed compatible with 
the pre-shock velocities estimated at the end of Section im 
for the Aab pair. If an important fraction of the X-ray over¬ 
luminosity was attributed to CWAB, the lack of significant 
variation between 2002 and 2014 would suggest the obser¬ 
vations actually occurred at somewhat symmetrical orbital 
phases, characterized by similar separations (and accord¬ 
ingly similar X-ray emission). This would however not be in 
agreement with the estimated orbital phases, and it is very 
unlikely that the uncertainty on the ephemeris of the wide 
orbit could lead to such a poor estimate of the orbital phases. 
Alternatively, this result may point to a lack of 1/D depen¬ 
dence of the the X-ray emission from CWAB, even though 
it would be at odd with the adiabatic nature of the shocked 
winds in the wide orbit. It is also interesting to note that 
the bright synchrotr on radio emission along the wide orbit 
jBlomme et al.ll2007l l provides compelling evidence that the 
shocks of the colliding-winds are clearly present. Otherwise, 
no particle acceleration could be operating. A wind-wind 
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interaction region does not therefore need to be bright in 
X-rays to be efficient at accelerating particles. This lends 
additional support to the idea that particle acceleration - 
and associated non-thermal radiative processes - can effi¬ 
ciently operate even in a wide range of conditions. 

Finally, now that approximate orbital elements are 
available, the marginal change of the X-ray flux between 
epochs A02 and S02 could certainly not be attributed to 
a change in the wide orbit. It could rather come from an 
eclipse-like event in the Aab close pair. On t he ba sis of the 
revised ephemeris given by llbanoglu et al.l il2013l ) , the or¬ 
bital phases in the Aab pair are 0.497, 0.922 and 0.446, re¬ 
spectively for epoch A02, S02 and S14. These numbers sug¬ 
gest indeed that the configuration of the system in Septem¬ 
ber 2002 was likely to lead to a partial eclipse by the primary 
stellar wind, although in April 2002 and September 2014 the 
system was closer to an eclipse by the secondary. However, 
considering the short period of the Aab pair and the time 
elapsed since the determination of its ephemeris, the uncer¬ 
tainty on the corresponding orbital phase may be too large 
to definitely ascertain this interpretation. 


5 SUMMARY AND CONCLUSIONS 

The X-ray emission from HD 168112 and HD 167971, two 
massive stellar systems in the open cluster NGC 6604, has 
been investigated on the basis of archive and new data ob¬ 
tained with the XMM-Newton satellite. Thanks to updated 
information about their multiplicity revealed in the past last 
years, the results of the observations could be discussed and 
interpreted in a more appropriate context. 

The X-ray spectrum of HD 168112 appears to be slightly 
overluminous, in agreement with the idea that a fraction 
of the emission comes from the colliding-wind region. It is 
moderately hard, with higher plasma temperature of 15- 
30 X 10® K. The analysis of the present X-ray time series 
allowed to estimate a conservative lower limit on the eccen¬ 
tricity of about 0.38. The latter estimate is based on the 
hypothesis of adiabatic shocked winds, which seems to be 
validated by a careful inspection of wind parameters along 
with recent estimates of the projected separation obtained 
using long baseline interferometry. 

The X-ray spectrum of HD 167971 is slightly softer than 
that of HD 168112. Thanks to the existence of an orbital so¬ 
lution, the X-ray spectral time series could be interpreted 
with a good idea of the orbital phases of the observations. 
Considering the ephemeris published on the wide orbit, the 
marginal variability measured between April and September 
2002 - if real - is more probably attributable to a partial 
wind eclipse in the short period Aab pair. The lack of sig¬ 
nificant variability between 2002 and 2014, despite the very 
different orbital phases, is certainly explained by a rather 
weak contribution of the X-ray emission coming from the 
colliding-wind region in the wide orbit (CWAB). The latter 
contribution is likely not bright enough to dominate the in¬ 
dividual stellar contributions and that coming from the close 
eclipsing binary (CWab). This is in agreement with the ra¬ 
diative regime that prevails for the shocked winds in the 
close binary, allowing for a high efficiency in the conversion 
of the kinetic power of the winds into thermal X-ray emis¬ 
sion. This efficiency is indeed expected to be significantly 


lower in the colliding-wind region in the wide orbit charac¬ 
terized by an adiabatic regime. This demonstrates that in 
hierarchical triple 0-type systems such as HD 167971, ther¬ 
mal X-rays do not necessarily constitute an efficient tracer 
to detect the long period wind-wind interaction. Interesting 
enough, a wind-wind interaction region does not need to be 
bright in X-rays to be efficient at accelerating particles. This 
lends additional support to the idea that particle accelera¬ 
tion can efficiently operate in a wide range of conditions. 

This study allowed to achieve a better view of the X-ray 
emission from these two systems. The present measurements 
will be of great importance when they will be confronted 
to theoretical models aiming at simulating the physics of 
colliding-winds in massive binaries. In particular, these mea¬ 
surements are important for energy budget considerations 
that will have to be addressed in the context of the inves¬ 
tigation of particle acceleration in massive binaries. Non- 
thermal processes - including particle acceleration and ra¬ 
diative mechanisms - occur on top the hydrodynamics of 
colliding-winds which is notably probed using X-ray mea¬ 
surements such as those discussed in the present study. 

From the observational point of view, campaigns are 
still in progress to derive the orbital elements of these sys¬ 
tems using the VLTI. On the other hand, hard X-ray investi¬ 
gations are envisaged to probe the non-thermal X-ray emis¬ 
sion from HD 168112 and HD 167971. As significant non- 
thermal radio emitters, the issue of the detection of non- 
thermal X-rays is still open. In this context, the information 
compiled in this study is intended to be used to prepare 
an overall X-ray investigation, including a soft thermal part 
(initiated in this study) and a hard non-thermal contribu¬ 
tion likely to be probed with NuStar or even ASTRO-H 
(expected to be launched in late 2015). 
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